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Abstract: We introduce a new approach to analyze the phase diagram of QCD at finite 
chemical potential and temperature, test it in the Gross-Neveu model at finite baryon 
density, and apply it to the study of the chemical potential-temperature phase diagram of 
QCD with four degenerate flavors of Kogut-Susskind type. 
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1. Introduction 

Hadronic matter at high density and finite temperature is expected to undergo one or more 
phase transitions from a confined phase to a quark-gluon plasma phase, and eventually a 
color superconducting phase due to a (weak) attractive channel between quarks of different 
colors, that could lead to the formation of quark pairs, analogous of Cooper pairing in solid 
state systems at low temperature. The behavior of hadronic matter at high density is also 
relevant in astrophysics, for the study of neutron stars and extremely compact objects, and 
in cosmology, for the understanding of phase transitions in the early Universe. Furthermore 
the new states of matter should be partially accessible to direct experimental verification 
in the heavy ion collision experiments. 

These features of finite density QCD can be studied only nonperturbatively. Unfortu- 
nately the lattice approach, the most powerful tool to perform first principles, nonpertur- 
bative studies, is affected in the case of finite density QCD by the well known sign problem, 
that has prevented for a long time any step toward its understanding. Indeed the moderate 
optimism generated by the finding of a first order phase transition at high density in the 
strong-coupling limit Q using the Glasgow reweighting method, was actually ill-founded 
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From the accumulated experience it was clear that the only way to minimize the effect 
of the phase of the fermion determinant in the integration measure of finite density QCD 
would be to work in the limit of large fermion masses. In fact taking the limit of fat QCD 
and using the MFA approach || in the numerical simulations, we were able to predict 
a tentative phase diagram for fat QCD jjj] which qualitatively confirmed the theoretical 
expectations. 

This was the rather unpleasant situation of three-color QCD at finite density until 
two years ago, when Philippe de Forcrand and Owe Philipsen ||, and later on D'Elia and 
Lombardo ||], reported some reliable results for light quarks, obtained from simulations at 
imaginary chemical potential where the sign problem, as well known, is absent. Fodor and 
Katz have also reported interesting results obtained within a modified version in two- 
parameter space of the Glasgow re-weighting technique |IJ] which, due to the re-weighting 
procedure, can be applied only to relatively small lattice sizes. A modified version of 
the last approach, which makes use of a Taylor expansion in the chemical potential, has 
produced also interesting results for larger lattice sizes in the small fi/T region [||. 

This paper is devoted to the presentation of a new approach to simulate QCD at finite 
temperature and baryon density, and to the application of this approach to the four-fermion 
model in three dimensions and to QCD with four degenerate quarks. Even if our theory 
scheme can resemble in some aspects the imaginary chemical potential approach, its range 
of applicability is at first sight much wider, as will become apparent after reading section 
2, which we devote to a detailed exposition of the approach. Section 3 contains the results 
of a test worked on the three-dimensional Gross-Neveu model at finite baryon density and 
in the large number of flavors limit. This model shows a non trivial phase diagram in 
the T, \x plane || and in addition, in the large N limit, can be analytically solved, thus 
allowing to confront results obtained within our approach with exact ones. In section 4 
we show our results for the phase diagram of QCD at finite chemical potential with four 
degenerate flavors of Kogut-Susskind type. The data reported in this section were obtained 
from simulations at values of the gauge coupling and fermion masses in the scaling region. 
Last section is devoted to collect our main conclusions and to discuss possible issues for 
near future. 

2. Theory scheme 

The lattice action for QCD at finite chemical potential fi and staggered fermions is 



S = S-p G + - ^2 i) n rji(n) (lln^n+i ~ ul_ iA ip n -i^j 

n i=l 

+ \ ^2 ^nVo{n) (e^U nfi ipn+o - e~ M [/^_ ^ n _ ) + m ^ tp n ^ n (2.1) 



2 

n 



where Spg is the standard Wilson action for the gluonic fields and rj^{n) the Kogut-Susskind 
phases. The temporal part of the fermionic action can also be written as 
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S T = COsh^-^VV^oW \U n fl1pn+0 ~ ^n-O.oV'n-oJ 
n 

+ sinh(^)i ^ it)o(ll) (Unfllpn+o + t^-o.o^n-o) (2.2) 

n 

The contribution proportional to cosh(/i) in 5 r is a fermionic bilinear denned through 
a matrix iA T with A T a hermitian matrix which anticommutes with 75. The contribution 
proportional to sinh(//) is also a fermionic bilinear defined through the hermitian matrix 
A T . Since the space-like part of the fermionic action is also a bilinear defined by the matrix 
ml + iA g , with A s hermitian, it is evident that the introduction of a chemical potential 
[i drive us to a complex fermion effective action. At the same time it is also obvious 
that an imaginary chemical potential solves the sign problem since cosh(/x) becomes the 
trigonometric cosine and sinh(/i) the trigonometric sine, but the last picking up a factor of 
i. 

2.1 The action 

Let us now define the following generalized action 

1 3 

S = S PG + -^2^4> n r]i( n ) (Un,iipn+i ~ Ul-i^n-i) + S T (x,y), (2.3) 

n i=l 

with 

S T (x,y) = x^^2ip n ri (n) (u n ,oi>n+o - L^_ 00 ^„_ ) 

n 

+ y\ ^2 ^o(n) (Unfl^Pn+O + f^-coVVi-o) (2.4) 

n 

where x and y are two independent parameters. The QCD action is recovered from S T (x, y) 
by choosing x = cosh(/x) and y = sinh(/i). 

Let us consider the generalized QCD action fl2.3j ). A natural question is: what should 
we expect for the phase diagram of this model in the x, y plane? 

We are going to conjecture a minimal phase diagram for this generalized model of 
QCD. To this end let us analyze first the y = line. The point x = 1 in this line would 
correspond to standard QCD at vanishing chemical potential. Now let be L t the temporal 
size of our lattice and assume we are in the scaling region but at a physical temperature 
T lower than the deconfining critical temperature T c . The point x = 1 in the y = 
line will be in the confined phase. If we increase now the inverse gauge coupling f3, the 
physical temperature increases and for (3 values large enough the point x = 1 will be in the 
unconfined phase. This strongly suggests the presence of a phase transition point in the 
y = line approaching the x = 1 point in this line by increasing (3 and eventually crossing 
it for values large enough. 
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Figure 1: Minimal phase diagram conjectured for the generalized model of QCD. The solid line is 
a line of phase transitions. The discontinuous line is the physical line x 2 — y 2 = 1. 

This argument can be enforced if we consider the strong coupling limit. In this limit, 
where the pure gauge action vanishes, one can change the temperature by using symmetric 
lattices and changing the temporal coupling in the fermion action. The point x = 1 in 
the y = line would be zero lattice temperature whereas x > 1 would correspond to 
a temporal lattice spacing a T smaller than the lattice spacing a in the space-like axis. 
Therefore x > 1 corresponds, in the strong coupling limit, to finite lattice temperature. 
But even in the strong coupling limit one reaches a chiral phase transition by increasing 
the lattice temperature flO|| , thus suggesting that at (3 = there is a phase transition point 
in the y = line for a critical value of x, x c , larger than 1. 

This argumentation drives one in a natural way to conjecture the minimal phase dia- 
gram shown in Fig. 1. The solid line would be a line of phase transitions that crosses the 
y = axis at values of x larger than 1 for small (3 values. By increasing (3 and keeping fixed 
the temporal lattice extent Lt, the critical point on the y = axis moves toward x = 1 
and eventually crosses it. The discontinuous line in this figure stands for the physical line 
x 2 — y 2 = 1 along which one recovers standard QCD. Along this line points with a coordi- 
nate y ^ would correspond to non vanishing chemical potential whereas x = 1 , y = is 
fj, = 0. The intersection of the solid line with the discontinuous one will therefore give us 
the critical chemical potential of QCD at a given temperature. A change in the physical 
temperature can be simulated by changing (3 keeping fixed Lt or vice-versa. In both cases 
the solid line in Fig.l will move and the intersection point which gives the critical chemical 
potential will change with the physical temperature. The nature of the phase transition 
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line will depend also on the other free parameters of the model i.e., number of flavors and 
quark masses. 

2.2 Symmetries of the partition function 

In order to get information on the symmetries of the partition function of the system 
described by the action ( |2,3| ), let us define the following change of variables in the path 
integral for the Grassmann fields, 

^j k = e ika ^ jk , jj k = e- ika ip jk , 

where k is the temporal coordinate of a given site, j stands for the three spatial coordinates, 
e ia is an element of the group Z^ t and Lt is the number of time slices. The only piece of 
the action Q2.3| ) not invariant under the previous variable change is S T (x,y), the different 
contributions of which transforming in the following way 



(^n%(«)^i_ ,o^n-o) = e~ ia ^ n r? o (n)f7^_ oo '0n-o) (2.5) 

In addition to the previous variable change for the Grassmann fields, we can also do a 
change of variables for the gluonic fields by multiplying each temporal link by an element 
e* 7 of the center Z% of the SU(3) group. Again the only non invariant contribution to (|2.3| ) 
is S T (x,y) and its different pieces transform in the same way as before, replacing e ia in 

Using these transformations we can conclude that the integrated partition function of 
generalized QCD has the following symmetry 

Z v (x + y, x - y) = Z v (>(x + y), e~ l \x - yfj , (2.6) 

where e ld is now an element of Z%L t . Since the partition function must be independent of 
the parameters introduced by changing variables in the path integral, the only possibility 
is that it depends on x, y through the combinations (x + y)(x — y), (x + y) 3fcit ; (x — y) 3fci t ; 
with k any natural number. Taking also into account CPT transformations we get that 
the dependence in the last two variables should appear as (x + y) 3kLt + (x — y) 3kLt , a 
result which could be also obtained by making the Polyakov loop expansion of the fermion 
determinant. Simple algebra tell us also that (x + y) 3kLt + (x — y) 3fcL * can be written as 
a combination of powers of (x + y) 3Lt + (x — y) 3Lt and of x 2 — y 2 . Therefore we can say 
without loss of generality that the partition function Zy(x,y) will be, at finite volume, an 
analytic function Zy(u,v) of the two variables 

2 2 

u = x — y , 

v = (x + y) 3Lt + (x- y f Lt , (2.7) 

and this is a non trivial information that we are incorporating in the description of the 
model. 
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Figure 2: Conjectured phase diagram in the x,y plane. We have incorporated the property of 
periodicity. The dashed line contains the only points accessible to numerical simulations of QCD 
at imaginary chemical potential. 



2.3 The extension to imaginary y 

The numerical analysis of the model ( |2.3[ ) for real values of x, y is not possible since, as 
stated before, we find again the sign problem. However if we take the y parameter as a pure 
imaginary number, y = iy, y real, numerical simulations are feasible since the fermionic 
matrix, as for the case of imaginary chemical potential, is the sum of a constant diagonal 
matrix plus an antihermitian matrix, the last anticommuting with the 75 Dirac matrix. 
The eigenvalues of the fermion matrix appear then as pairs of complex conjugate numbers, 
and the fermion determinant is therefore real and positive. 

Let be z = x + iy the complex number which we can associate to each pair of real 
values for x, y. We can write z as z = pe iri with 

p=(x* + f) 1/2 , tan 77=^. (2.8) 

x 

The previous analysis on the symmetries of the partition function tell us that for 
imaginary values of the y parameter, this function will depend only on p and p 3Lt cos(3L^) 
i.e., the free energy will be a periodic function of rj with period equal to 2ir/(3Lt). In 
particular if the phase transition line of Fig. 1 continues to imaginary values of y, the 
conjectured minimal phase diagram in the x,y plane would be that of Fig. 2, where we 
have incorporated the property of periodicity. 
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Figure 3: Phase diagram in the [u,v) plane. The solid line v = 2u? Lt corresponds to the y = 
line of Fig. 1. The physical line u — 1 and the hypothetical phase transition line are also included. 

In Fig. 2 we have also included the line p = 1 (dotted line) which contains the only 
points accessible to numerical simulations of QCD at imaginary chemical potential. As 
stated in the introduction, one can see now how this approach has in principle much more 
potentialities than the imaginary chemical potential approach. Indeed by increasing the 
inverse gauge coupling (3, the phase transition line of Fig. 2 moves approaching more and 
more the origin of coordinates. In some interval (/3 m , Pm) the transition line intersects the 
p = 1 line and then a phase transition will appear at imaginary chemical potential. In 
such a situation, the physical temperature is so high that the system is in an unconfined 
phase for any real value of the chemical potential. In order to understand something on the 
actual phase diagram from simulations at imaginary chemical potential, we need not only 
to extrapolate from imaginary to real values of p but also to extrapolate from higher to 
lower temperatures or equivalently from higher to lower values of (3. The advantage of our 
approach is that we can avoid the second kind of extrapolations because in our simulations 
p is not enforced to be 1. 

This discussion can become more clever if we go from the x, y plane to the u, v plane, 
where u and v are the variables denned in ( |2.7| ), The reason for choosing this plane is that 
it contains simultaneously the two planes of Figures 1 and 2. In Fig. 3 we have plotted 
the line 

v = 2vh u (2.9) 

which corresponds to the y = line of Fig. 1. We have also included in Fig. 3 the line u = 1 
which contains the discontinuous or physical line of Fig. 1 (points up the intersection point 
with the line v = 2u^ Lt ) as well as the dotted line of Fig. 2 (points down the intersection 
point). All the plane of Fig. 2 maps into the region down the line v = 2u^ Lt of Fig. 3, 
and this region is therefore accessible by simulating at imaginary y. On the other hand, 
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the plane of Fig. 1 maps into the region up the solid line of Fig. 3. The only points in 
this Figure accessible to numerical simulations at imaginary chemical potential are those 
on the line u = 1, down the intersection point with the line v = 2u^ Lt . 

In Fig. 3 we include also a hypothetical but plausible realization of a phase transition 
line at a given physical temperature (simulations would be performed at fixed (3 and Lt). 

2.4 Determination of the critical chemical potential near the continuum 

In order to get the critical value of the chemical potential we should determine the coordi- 
nates of the intersection point between the solid line and the physical line in Fig. 1. From 
the symmetries of the partition function we know that the phase transition line is an even 
function of x and y. We can write therefore the following equation for the critical line of 
Fig. 1. 



1 + a (ft L t ) + a 2 (ft L t ) y 2 + a 4 (ft L t ) y A + 0(y e 



(2.10) 



By fixing the lattice temporal extent Lt and the gauge coupling (3 one fixes the physical 
temperature T. The intersection point of the critical line with the physical line 



y 2 = a (ft L t ) + a 2 (ft L t ) y 2 c + a 4 (ft L t ) y\ + 0(y 6 c ), (2.11) 

will give us the critical value of the chemical potential (y c = sinh(/x c o)) at this temperature. 

If we approach the continuum limit by fixing the physical temperature T, i.e, by sending 
(3 to oo keeping the product of L t times the lattice spacing a fixed, we can write equation 
( |2.11 ) in the following way 



1 it „ 
+ 3#° + 



a (ft L t ) L\ + a 2 (ft L t 



+ a 4 (ft L t ) + -a 2 ((3, L t ) 



L: 



+ ... 



(2.12) 



where in order to get this equation we have multiplied both sides of (2.11) by L 2 . The 
coefficients (ft Lt) will depend in general on Lt and on the dimensionless parameter A/T, 
where A is the dimensional renormalization group parameter that fixes the QCD scale, and 
we have neglected for simplicity in all this analysis the dependence on the quark masses. 
Equation ( |2.12 ) then can be rewritten as follows 



- "n ( ^,L t ) Lf + I a> 



+ 



A \ 1 /A r 
CL4, | — ,L t I + -02 ( — ,Lt 



1 



a 2 L 2 



L\ 



-2 / M< 



T J 3 V T 7 3 " 
which suggests the following scaling laws for the coefficients 



+ 



(2.13) 



A _ 



r, \ - J <(„ \-,L t ) L k t , 



(2.14) 
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The critical value of fi/T will be given by the solution of the equation ( |2.13| ), This 
equation suggests that for fJ, c /T c small, the determination of the first two coefficients of 



(2.10) could be enough for the computation of the critical chemical potential. The goodness 



of this approximation will depend on the values of the coefficients of (2.10). In the next 



section we show that this approximation works quite well in the Gross-Neveu model, even 
at large values of fi c /T c . 

The strategy for the determination of the critical chemical potential is then the follow- 
ing. From numerical simulations at imaginary values of y, y = iy, near the phase transition 
point ((l + a ) 1/2 ,0) one can locate several phase transition points in the (x, y) plane (see 
Fig. 2). By fitting these points with equation ( p.lOf ) with the + sign of the coefficient 
proportional to y 2 replaced by -, we can numerically measure the first coefficients. The 
critical value of the chemical potential, fj, c , will then be given by 

\ 1/2 

a x 



Mc = ±sinh-M I . (2.15) 

An alternative approach for the determination of the critical chemical potential is to 
use the variables (u,v) instead of (x,y). Since near the continuum u takes values near 1, 
the phase transition line of Fig. 3 could be fitted with the following equation 

v = b + h (u - 1) + b 2 (u - l) 2 + O ((« - l) 3 ) . (2.16) 
The coefficient bo fully determines the critical [i through the following equation 

^ = _L cosh -i ( h\ (2.17) 



3L t V 2 , 

The use of the (u, v) variables incorporates information on the symmetries of the free 
energy density. Furthermore consistency between results obtained from fits with equations 
(2.10) and ( 2.16j ) would improve the confidence level on the full procedure. 



2.5 Two other alternative ways for extracting the critical /x: the (f3,y) and ((3,x) 
planes 

In the previous paragraph it was assumed that the inverse gauge coupling f3, the quark 
masses and the lattice temporal size Lt were fixed and then we explored the phase diagram 
in the (x, y) plane. In physical terms this means that, if (5 is large enough to be in the 
scaling window, we fix the physical temperature T and the physical quark masses. However 
it could be interesting in practice to analyze the projection of the phase diagram into the 
other two planes {(5, y) and (/?, x) for the following two reasons: 

i. it would give an independent determination of (/u c ,T c ); 

ii. one could apply some powerful techniques in the numerical analysis, as the Ferrenberg- 
Swendsen approach [11]. 
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beta 



Figure 4: Conjectured phase diagram in the ((3, y) plane. The solid line is the physical line whereas 
the discontinuous one is a line of phase transitions. 

Let us start with the case in which we fix x = xo > 1, Lt and the lattice quark masses. 
The conjectured phase diagram in the (/?, y) plane is reported in Fig. 4. The solid line is 
the physical line 

y = yo = sinh (cosh _1 (x )) , 

and the discontinuous line is a line of phase transitions. Using again the symmetries of the 
partition function we can write for this line the following equation 

P = po (x ,Lt) + b 2 (xo, L t ) y 2 + (y 4 ) , (2.18) 

where we have neglected for simplicity the dependence of the coefficients on the quark 
masses. 

We are interested in the determination of the intersection point between the physical 
and phase transition lines of Fig. 4. As in the case discussed in the previous section, one 
can measure (3q and 62 from simulations at y = and at imaginary y = iy (keeping xq 
fixed). Once we estimate the lattice spacing a, assumed we are in the scaling region, we 
can obtain /3 C , \i c and the critical values of the physical quark masses. Notice also that in 
the imaginary chemical potential approach xo is enforced to be less than 1. 

The last case to analyze, the projection of the phase diagram into the (j3, x) plane, is 
less interesting from a practical point of view because: 

i. The physically relevant region in the (/3, x) plane (keeping y = yo fixed) is x larger but 
near to 1. The equation for the phase transition line near x = 1 can therefore be written 

as 

P = Po (Vo) + ai (y ) (*-!) + a 2 (y ) (x - l) 2 + O ((x - l) 3 ) . (2.19) 
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The critical inverse gauge coupling (3 C will be obtained by putting x = xq in ( 2.19 ), with 



xq = cosh (sinh (yo)) ■ 

But since no symmetry enforces now a\(yo) to vanish, we should do, in contrast with the 
previous cases, a three parameters fit to get j3 c to order (xo — l) 2 . 

ii. Numerical simulations must be performed at imaginary values of y and therefore we 
need, at the end of the procedure, to do also an extrapolation of (3 c {y) from imaginary 
values of y to real y = yo- 

2.6 Vacuum expectation values on the phase transition line 

In order to get a deep understanding of the phase transition(s) at finite baryon density 
we would like to know not only the value of the critical chemical potential at a given 
physical temperature but also how vacuum expectation values of physical quantities change 
when crossing the phase transition line. To this end we will focus here on two typical 
observables, the mean plaquette energy and chiral condensate, and on the first theory 
scheme described previously i.e., the analysis of the phase diagram in the x,y plane at 
fixed (3. The generalization to other observables as the baryon number density, and, or to 
other projections of the three-dimensional phase diagram will be straightforward. 

The chiral condensate q(x 2 ,y 2 ) and plaquette energy E(x 2 ,y 2 ) at fixed (3 will be an- 
alytical even functions of x and y in all the points very near to the phase transition line 
both on the left and on the right of it. If we compute for instance the chiral condensate 
q(x 2 ,y 2 ) in all the points of the phase transition line 

x 2 = 1 + a (A L t ) + a 2 (f3, L t ) y 2 + a 4 (/?, L t ) y A + 0(y 6 ) , 

approaching it from the left-hand side, we expect an analytical even function of y. Then 
we can write for the chiral condensate along this line the equation 

Qc- (y 2 ) = q?2 + £V + 0_ (y 4 ) , (2.20) 
and the same for all the points on the right of the phase transition line 

Qc + (y 2 ) = 5$ + & + + (/) . (2.21) 

From simulations at imaginary y = iy one can measure the chiral condensate in both 
sides along the phase transition line, and by fitting the numerical results with these equa- 
tions, we would get a numerical determination of the parameters in ( 2.20| ), (2.21). This 



procedure would allow us to extend the chiral condensate results to the physical point 

y = y c - 

3. The Gross-Neveu model at finite density 

To check all the conjectures and ideas developed in section 2, we decided to use as a toy 
model the four-Fermi model in three dimensions with discrete chiral symmetry at finite fi 
and T. The continuum Lagrangian is 
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L 



N f 

£ 

3=1 



2N 



f 



^O')^(i) 



(3.1) 



where i/j^ is a four-component spinor, Nf is the number of flavors and the discrete Z2 
symmetry is ip — > 75?^, t/> — > —'075. This model has a non trivial renormalization group 
fixed point characterized by the spontaneous breaking of the chiral symmetry at a finite 
critical coupling g^, its 1/Nf expansion about the fixed point g^ being renormalizable. In 
addition the model can be solved analytically in the large Nf limit and then we can compare 
numerical predictions with exact results in this limit. Finally the phase structure of the 
model in the (/i, T) plane resembles to that expected for QCD, with a phase transition line 
separating a phase where chiral symmetry is spontaneously broken from a phase where the 
symmetry is restored. 



After the introduction of an auxiliary scalar field a, equation (3.1) becomes 



Nf 



(3) 



3=1 



N f 2 



(3.2) 



the vacuum expectation value of a being a good order parameter for the chiral symmetry 
restoration transition. 

Regularized in a space-time lattice with staggered fermions, the Gross-Neveu action is 

H 



Nf/2 3 1 

S = E E 2 r] ^ n ^ i ^ ^ n + ~ 

i=l n f_i=l 

1 N,/2 N 

+ g E E^( n )^( n ) E + E a2 (")' ( 3 - 3 ) 

i=l n <n,n> n 

where the symbol < n,n > denotes the set of eight dual lattice sites h surrounding the 
lattice site n, rjo(n) = 1 and i]^(n) = (— l) ra f+ n M-i for fj, = 1, 2. 

The inclusion of a chemical potential /j, selecting a finite fermion density is achieved, as 
in QCD, by multiplying all forward (backward) temporal links by e^, (e - ^). The fermionic 
bilinear part of the action then becomes 

Nf/2 

s F = E 2^ (n) [ e ^ i(n + 0) " e ~^ n - °)1 

i=l n 
JV>/2 2 

+ E E E ^ii^Vain) [M n + a) - - a)] 

j=l n a=l 
JV>/2 

+ gEE*( n W ft ) E *(*o ( 3 - 4 ) 

i=l n <n,n> 
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Figure 5: Phase diagrams of the generalized four-Fermi model in the (x,y) (5a), (u,v) (5b) and 
(l/g 2 ,y) (5c) planes respectively. All the phase diagrams correspond to a lattice temperature 
1/L t — 1/4 in lattice units. The inverse four- Fermi coupling 1/g 2 = (3 = 0.75 in Figs. 5a and 5b 
whereas x — 1.0211 in Fig. 5c. The solid line is the physical line whereas the discontinuous line is 
a second order phase transition line. 
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In the large number of flavors limit, Nf — > oo, fluctuations around the saddle point 
solution are suppressed and then the vacuum expectation value £ of the scalar field a can 
be self-consistently determined using the gap equation 

S = -g 2 , (3.5) 

where Otpi/)} here means the mean value of the chiral condensate in the free fermion theory 
at finite chemical potential and with an effective fermion mass m = X. 

The formalism developed in the previous section can be straightforwardly applied here. 
Indeed we can define a generalized fermion action 



S F (x,y) = x- ^2 ^2iJi(n) [tpi(n + 0) - ^(n - 0)] 

i=l n 
l N f /2 



2 

i=l n 



N f /2 2 

+ ^2 ^2 ^2 nM n )ria(n) [ifji(n + a)-ipi(n- a)] 



2 

i=l n a=l 



Nf/2 



j=l n <n,n> 



that agrees with action ( |3.4| ) for x = cosh(/i), y = sinh(/x). In the large iVj limit the 
generalized model can also be analytically solved, and the solution is again given by the 
gap equation ( |3.5| ) where now (ipip) means the mean value of the chiral condensate in the 
generalized free- fermion model at an effective fermion mass m = E. 

We are now in good conditions to check if the conjectures formulated in section 2 are 
realized in the four-fermion model. Figures 5a, 5b and 5c show the phase diagrams of the 
generalized four-Fermi model in the (x,y), (u,v) and (y,g 2 ) planes respectively. All the 
phase diagrams correspond to a lattice temperature 1/Lt = 1/4 in lattice units. The inverse 
four-Fermi coupling j3 = 1/g 2 = 0.75 in Figs. 5a and 5b whereas x = 1.0211 (fi = 0.2) 
in Fig. 5c. A second order phase transition line separating a broken from an unbroken 
phase and crossing the physical line x = cosh fj,, y = sinh \x was found in all the cases, thus 
confirming actually all our conjectures. 

To have a guess on the potentialities of our approach when applied to QCD we have 
done a further check. This check consisted in fitting several phase transition points in 
the (x,y) plane at imaginary y with equation ( 2.10| ) in the quadratic approximation (two- 



parameter fits). The purpose of the check was to assume these phase transition points as 
they were obtained from numerical simulations at imaginary y, where the sign problem is 
absent, and then to check the degree of precision with which we can predict \x c and T c . 

To locate the physical region in which the agreement between the extrapolated results 
and the exact ones is better, we have plotted in Fig. 6 the phase diagram in the (n/T, T/So) 
plane for several values of the lattice temporal extent Lf. Eq is the expectation value of 
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Figure 6: Phase diagram in the (/z/T, T/£o) plane for several values of the lattice temporal extent 
Lt- So is the expectation value of the scalar field a in the ground state at T = 0. The lines 
correspond to analytical results at the same values of Lt while the symbols correspond to the 
extrapolated results. 



the scalar field a in the ground state at T = 0, and the plots have been done assuming 
asymptotic scaling. The lines correspond to analytical results. The observed dependence 
on Lt, especially at large fJ,/T, shows that the effects of regularization are still important 
in this region for these lattices. One can also observe in this figure how the extrapolated 
results agree very well with the exact ones even at large values of n/T. 

Finally and in order to confront our method with the imaginary chemical potential 
approach, we have plotted in Fig. 7 the phase diagram of the Gross-Neveu model in 
the fi,l/g 2 plane as obtained from the analysis of the x, y plane at imaginary y in the 
quadratic approximation of equation ( 2.1C| ) (crosses), and from the imaginary chemical 



potential approach (dashed lines), the last obtained from quadratic as well as quartic fits 
of the critical coupling as a function of the imaginary chemical potential. The solid line 
stands for the exact results on a lattice of temporal extent Lt = 4. As can be seen both 
approaches agree with the analytical results in the small chemical potential region, but 
results obtained with our method follow quite well the exact ones in the large \i region in 
contrast with the imaginary chemical potential results that give even negative values for 
g 2 at large \i in the quadratic approximation. From both theoretical and practical points 
of view there is little hope to get reliable values for more than two-three coefficients from 
imaginary \i numerical data, so we improved the fit procedure, as stated before, considering 
up to quartic fits of the critical coupling as a function of the imaginary chemical potential. 
As can be seen in Fig. 7 the extra degree of freedom improves the results at small \x, but 
not in a relevant way. 
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Figure 7: Phase diagram of the Gross-Neveu model in the plane as obtained from the 

analysis of the x,y plane at imaginary y (crosses), and from the imaginary chemical potential 
approach (dashed line). The solid line stands for the exact result on a lattice of temporal extent 
L t = A 

4. Four flavor QCD at finite T and fx 

Understanding the finite temperature and finite density phase structure of QCD is one of 
the main challenges of lattice field theory. Whereas qualitative and quantitative progress 
on the phase structure of the model at vanishing baryon density has been achieved in the 
last years, the sign problem has prevented, as discussed in the introduction of this paper, 
any significant progress on the understanding of the behavior of matter at high density. 

We want to devote this section to the application of the approach described in the 
previous sections to the analysis of the phase diagram of four flavor QCD at finite baryon 
density. The aim of such an analysis is twofold: first we want to verify if the conjectures 
formulated in the first section, and corroborated in the previous section for the Gross-Neveu 
model, do actually work in QCD, and second we want to gain some experience in order to 
analyze in a future work the more interesting cases of 2 or 2+1 flavor QCD. 

Even if a good understanding of the phase diagram of the four flavor model would be of 
interest, the aim of this section is not to do an exhaustive analysis of the four flavor model 
but rather to verify if the approach reported in this paper can hopefully be applied to QCD. 
This is the reason why we have concentrated our efforts on the analysis of one point of the 
phase diagram, but using two of the independent approaches described in section 2: the 
determination of the critical values of fi and T from the analysis of the phase diagrams in 
the (x, y) and {(5, y) planes. Since the two approaches are based on independent numerical 
simulations and independent fits, the agreement between the predicted critical parameters 
should be considered a proof of the reliability of the results. 
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4.1 The (x,y) plane 

The model analyzed is therefore that described by action ( |2.1[ ) i.e., QCD at finite \x and 
four degenerate staggered flavors. In order to determine a point of the phase transition 
line in the (T, [i) plane we performed HMC simulations of the generalized model described 
by action (|2.3[), ( |2.4[ ) at several values of x and imaginary values of y = iy. The lattice 
size was 8 3 x 4, the inverse gauge coupling (3 = 5.05, and the quark masses were fixed to 
ma = 0.06 in lattice units. The corresponding physical temperature is T ~ 136MeV. 

We have performed simulations for 7 values of y from to 0.15; for each value of y 
8 — 13 values of x have been simulated with statistics ranging from 10000 to 50000 HMC 
trajectories (the larger ensembles near the critical points). The simulations have been 
performed on the Linux clusters of LNGS-INFN, using a total of 50 CPUs. 




Figure 8: Measured critical points in the (y,x) plane at /? = 5.05 and ma — 0.06. The upper and 
lower half figure correspond to imaginary and real y respectively. The solid lines are a quadratic fit 
of the points and its analytical continuation to real y. We include also the physical line x 2 — y 2 = 1. 
Dashed lines represent maximum deviation of the fitting function within statistical errors. 

From these numerical simulations we have determined seven critical points in the (y, x) 
plane (see Fig. 8), looking for the position of the maximum of the plaquette, chiral, and 
Polyakov loop susceptibilities. The position of the maximum was compatible for the three 
observables in all the cases, but the results for the plaquette susceptibility seem to be of 
better quality. The two-state signal observed for the three measured quantities in all the 
points plotted in Fig. 8 indicates a first order phase transition line. Unfortunately we have 
no data on larger lattices which would allow a finite size scaling analysis to confirm these 
expectations. 
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By fitting the seven points of Fig. 8 with a quadratic function (equation ( 2.10| ) with 
a n = for n > 2) we get the solid line in this Figure. Its analytical continuation to real 
values of y has also been plotted in Fig. 8 where one can also observe the intersection 
point with the physical line x 2 — y 2 = 1. The coordinates of this intersection point x c = 
1.0081(2), y c = 0.127(2) allow us to predict that the critical chemical potential in lattice 
units is fi c = 0.127(2) at j3 = 5.05 and ma = 0.06. 

We can fix the lattice scale from the mass of the p resonance, 770 MeV. In such case the 
results of reference [12] imply that the lattice spacing is approximately 0.3 fm at /3 = 5.2 
for ma = 0.06, and we set the lattice spacing at smaller values of (3 by using the two 
loop j3 function. The phase transition at \i = is located at f3 ~ 5.08 [O]. This gives 
Tq ~ 136 MeV. Our results reported in the previous paragraph give another point of the 
phase transition line: Tq 133 MeV, fic ~ 68(1) MeV. 




Figure 9: Measured critical points in the (y, 0) plane at x = 1.009 and ma = 0.06. The upper and 
lower half figure correspond to imaginary and real y respectively. The solid lines are a quadratic 
fit of the points and its analytical continuation to real y. We include also the physical line y = 
sinh(cosh _1 (1.009)). Dashed lines represent maximum deviation of the fitting function within 
statistical errors. 



Once we have determined the coordinates of the critical point, we can apply the for- 
malism described in section 2.6 to analyze the nature of the physical phase transition. A 
very interesting feature that emerges from our simulations is that the latent heat is inde- 
pendent of the parameters x,y within errors in such a way that there is not a real need 
to perform an analytical continuation; therefore the physical phase transition is also first 
order. Similar results have also been obtained for the chiral condensate and Polyakov loop. 
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4.2 The (f3,y) plane 

By fixing x = 1.009 (corresponding to \i = 0.134, our first rough estimation of the critical 
value at (3 = 5.05 and ma = 0.06), we have also performed simulations in the (/3,y) plane, 
using statistics and number of values for the parameters (3 and y similar to those of the 
(x, y) case. The aim of this calculation, as stated before, is to check if we get results 
consistent with those of section 4.1. 

Fig. 9 is the same as Fig. 8 but in the (y, j3) instead of the (y, x) plane. The solid lines 
are a fit of the measured phase transition points with equation (|2.18[) and its analytical 
continuation to real values of y. Again in this case a clear two state signal was observed for 
the three measured observables. As in Fig. 8, we have also plotted in Fig. 9 the physical 
line that in this case is the line y = sinh(cosh _1 (1.009)). The intersection point of the 
phase transition line with the physical line gives the critical value j3 c = 5.046(1) in close 
agreement with the results reported in the previous section, if we take in mind that the 
value of \i used (0.134) is slightly larger than the actual critical (j, at = 5.05 (0.127). For 
this case the critical point is Tq sa 132 MeV, (jlq sa 71 MeV. Also for this set of simulations 
the global behavior of latent heat suggests a first order transition. 

5. Conclusions and comments 

We have introduced a new approach to analyze the phase diagram of QCD at finite chemical 
potential and temperature, and applied it to the four-fermion model in three dimensions 
and to QCD with four degenerate quarks. 

Our scheme is based on the introduction of a generalized QCD action which depends 
on two extra free parameters, x, y. When the point (x, y) belongs to the curve x 2 — y 2 = 1 
one recovers standard QCD at finite chemical potential. 

The generalized QCD action suffers also, for real values of y, from the sign problem. 
However at imaginary values of y the fermion determinant is positive definite and therefore 
standard numerical simulations are feasible. 

Even if our theory scheme can resemble the imaginary chemical potential approach, 
we have shown that its range of applicability is much wider. Indeed our approach shares 
with the imaginary chemical potential method the need of analytical. However whereas 
imaginary chemical potential simulations have to be performed at high temperature (large 
(3) i.e., in the deconfined phase for real values of fj,, in our scheme simulations can be 
performed at any temperature and any value of fi/T. The price to pay for that is one more 
parameter in the numerical procedure. 

To make evident the potentialities of our approach we have confronted it with the 
imaginary chemical potential method when applied to the Gross-Neveu model in the limit 
of large number of flavors, where analytical results are available. We have shown how both 
approaches agree with the analytical results in the small chemical potential region, but 
results obtained with our method follow quite well the exact ones in the large fi region in 
contrast with the imaginary chemical potential results, that give even negative values for 
g 2 at large fj,. An open point however is to dispose of a good criterion to understand in 
which region of the parameters space the quadratic approximation of equation ( 2.1C| ) (or 
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equivalently of equation ( 2.13 )) works well. For what concerns the Gross- Neveu model, the 
aforementioned region is actually very wide. 

We have also performed simulations of QCD with four degenerate flavors in order to 
test the procedure in a real QCD simulation. We have determined the critical chemical 
potential for a value of (3 inside the scaling region and, using the second scheme depicted 
in section 2.5, we found the critical f3 for a given /i value, obtaining a consistent pair of 
results. 

As stated before, our theory scheme requires more computer time than other more 
standard simulations since we have one more parameter in the action. Roughly speaking 
the computer time needed to get a single point of the phase diagram in the (/i, T) plane 
is of the order of the one needed to get the critical line within the imaginary chemical 
potential approach. However, our approach, as we have shown for the Gross-Neveu model, 
seems able to reproduce quite well the phase diagram at values of fj,/T not accesible to the 
other standard methods, and this is in our opinion the reason that makes it relevant for 
the field despite of its computing cost. 

Concerning the Glasgow reweighting procedure in two parameters space , one could 
consider the possibility of using the fermion determinant at the critical point of the y = 
line instead of the determinant at ^ = in the integration reweighted measure. We suspect 
that hereby one could improve the overlaps. 

Finally we plan for near future to apply our approach to the analysis of two-color QCD 
at finite \x and to 2 and 2+1 flavor QCD. 
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